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HIGHLIGHTS 


•  Ceramics  were  prepared  by  conventional  and  uniaxial  pressing. 

•  High  vapour  pressure  was  demonstrated  to  achieve  good  proton  conduction. 

•  The  solid  electrolyte  obtained  showed  an  excellent  conductivity  below  600  °C. 

•  Modelling  of  E1S  spectra  allowed  separating  bulk  and  grain-boundary  conductivities. 

•  Proton  and  anionic  conduction  through  the  solid  electrolyte  were  separated. 
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Because  of  their  high  ionic  conduction  even  at  relatively  low  temperatures,  proton  conducting  ceramics 
are  one  of  the  most  promising  electrolytes  for  fuel  cell.  In  contrast  to  oxide-ion  conductors,  proton¬ 
conducting  systems,  especially  in  an  electrolyser  plant,  could  operate  below  600  °C,  critical  tempera¬ 
ture  for  mechanical  and  hot  corrosion  resistance  of  common  stainless  steels. 

Proton  conduction  in  perovskite  type  ceramic  was  analysed  under  the  water  molecules  insertion  to 
maximize  proton  conduction.  To  this  end,  the  SrZro.9Lno.1O2.95  electrolyte  and  its  microstructure, 
particularly  grain  size  was  investigated.  The  insertion  and  release  of  water  in  the  perovskite  structure 
was  followed  by  impedance  spectroscopy  and  by  using  the  so-called  Brick  Layer  Model  (BLM)  to  separate 
the  bulk  and  grain  boundary  conductivities.  It  was  found  that  the  conductivity  in  both  zones  increases 
namely  with  steam  pressure.  It  was  also  verified  that  the  bulk  conductivity  is  higher  than  that  of  grain 
boundaries.  Consequently  both  the  operating  steam  pressure  and  the  grain  size  appear  as  the  main 
parameters  that  can  be  tuned  to  enhance  the  proton  conduction. 
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1.  Introduction 

Hydrogen  could  play  an  increasingly  important  role  as  a  future 
energetic  carrier.  Its  demand  will  undergo  a  substantial  increase  in 
the  future  as  an  energy  resource  for  transport,  especially  for  fuel 
cells  or  for  synthetic  hydrocarbons  used  by  the  internal  combustion 
engine.  Other  needs  will  come  from  chemical  industry,  mainly  in 
terms  of  ammonia  production,  and  also  for  the  treatment  of  heavy 
hydrocarbon  resources  thereby  extending  the  life  of  oil  fields  [1—3], 
As  a  result,  mass  production  of  hydrogen  at  a  competitive  cost 
becomes  a  challenge,  but  this  production  must  avoid  greenhouse 
gas  formation.  This  contrast  to  the  steam  reforming  largely  used  in 


0.  Lacroix  et  al.  /  Journal  of  Power  Sources  270  (2014)  506—515 


507 


industrial  scale,  but  produces  CO2.  The  electrolysis  of  water  meets 
these  criteria;  especially,  high  or  intermediate  temperature  elec¬ 
trolysis  (HTE  or  ITE)  which  offers  economical  prospects  compared 
with  classical  electrolysis.  The  energy,  electricity  and  heat,  neces¬ 
sary  for  the  electrolysis  process  will  be  provided  by  nuclear  elec¬ 
trical  power  plants  [4], 

Currently,  the  water  electrolysis  at  high  temperature  in  liquid 
phase  operates  an  energy  efficiency  of  92%.  HTE  or  ITE  offers  the 
following  benefits: 


Up  to  now,  trivalent  cation  —  substituted  perovskite-type  oxides 
[9—11,16—29,41—47]  were  reported  to  show  rather  good  proton 
conduction  above  450  °C.  Based  on  these  results,  the  objective  of 
this  paper  was  to  clarify  the  role  of  water  insertion  into 
SrZro.gLno.i 02.95  [48]  substituted  by  lanthanide  element  (Ln)  to 
optimize  proton  conduction  with  respect  to  perovskite  composi¬ 
tions  and  sintering  conditions  [49—54], 

2.  Brick  layer  model 


-  Decrease  the  thermodynamic  equilibrium  potential  (£0)  of  wa¬ 
ter  decomposition  process, 

-  Decrease  the  activation  overpotential  at  the  electrode, 

-  Decrease  the  ohmic  drop  inside  the  electrolyte, 

-  Use  of  the  solid  electrolyte  needs  no  diaphragm,  avoiding 
corrosion  and  leaks. 

One  of  the  current  objectives  of  energy  production  industry  is 
the  development  of  an  electrolyser  able  to  provide  H2  in  the  in¬ 
termediate  temperature  range  400—600  °C,  to  ensure  a  good  safety 
for  the  material.  In  this  temperature  range,  common  stainless  steel 
can  be  used  for  the  structure  of  electrolysing  cell  instead  of 
expensive  nickel  rich  alloys. 

According  to  literature  [5—8],  satisfactory  ion  conductivities 
could  be  reached  with  proton  conducting  ceramics  such  as 
substituted  perovskites,  in  this  temperature  range.  Their  properties 
were  studied  for  the  first  time  by  Iwahara  and  his  group  [9—11], 
Since  their  pioneer  investigations,  many  works  on  perovskite  as 
proton  conducting  electrolyte  are  reported  in  recent  years  for  in¬ 
termediate  temperature  fuel  cell  or  steam  electrolysis  [12—29], 
The  perovskite  with  representative  formula  ABO3,  has  a  cubic 
structure.  The  B  sites  are  tetravalent  cations  (Zr,  Ce,  Nb,  Ti . . . ).  The  A 
sites  are  generally  divalent  cations  (Sr,  Ba,  Ca  ... )  which  form  planes 
between  BC>6  octahedrons.  These  atoms  have  high  ionization  ten¬ 
dency,  therefore  a  long-range  interaction  is  expected.  Proton 
insertion  occurs  for  oxygen-deficient  structures  exposed  to  high 
temperature  water.  Oxygen-vacancies  are  generally  formed  by 
substituting  the  stoichiometric  perovskite  ABO3  with  lower- 
valence  ions  D  to  the  B-site.  D  sites  are  trivalent  cations,  often 
lanthanide  Ln  (Ce,  Sm,  Ho,  Er,  Yb, ...).  As  a  whole,  the  ABi_*Lnx03_x 
compound  is  formed. 

This  species  is  not  intrinsically  proton  conductor;  therefore  it  is 
necessary  that  oxygen-deficient  perovskite  will  interact  with  water 
molecules  of  the  steam  phase.  It  is  assumed  that  the  latter  is 
dissociated  into  hydroxyl  ions,  filling  oxygen  vacancies  (Eq.  (1)). 

Og  +  Vo  +  H20^20Hq  (1) 

Therefore,  the  number  of  charge  carrier  sites  depends  on  both 
the  number  of  oxygen  vacancies  sites  and  on  the  amount  of  inser¬ 
ted  water.  The  mechanism  for  proton  transfer  in  these  oxides  re¬ 
mains  yet  unclear  due  to  the  difficulty  of  detecting  proton  motion 
directly  through  experiment.  It  is  recently  shown  that  the  adsorp¬ 
tion  of  some  species  like  (hydroxide,  hydrocarbonates)  avoid  the 
incorporation  of  bulk  protons  [30,31], 

The  conductivities  of  electrolytes  are  dependent  on  the  proton 
mobility  and  proton  concentration.  Some  correlations  can  be  found 
in  literature  between  electrochemical  properties  of  electrolytes  and 
the  nature  of  constituent  elements.  It  should  be  noted  that  not  only 
the  role  of  constituent  elements  of  the  ABO3  perovskite  [32,33],  but 
also  the  nature  and  mass  of  the  elements,  the  net  alkalinity,  the  size 
and  concentration  of  substituent  elements  [34,35]  on  the  perfor¬ 
mances  are  reported.  In  the  choice  of  an  electrolyte,  the  sensitivity 
to  carbon  dioxide  and  the  moisture  resistance  [36—38]  need  also  to 
be  considered  [39,40], 


For  modelling  the  conductivity  of  polycrystalline  ceramic  ma¬ 
terials,  a  brick  layer  model  (BLM)  generalised  by  Nafe  [55]  was  used 
[49,50],  and  illustrated  in  Fig.  1. 

In  this  model,  the  bulks  (grains)  are  assumed  to  be  of  cubic 
shape  of  edge  length  C  surrounded  by  a  thin  grain  boundary  layer 
of  thickness  g.  The  resistance  (R)  and  the  capacitance  (C)  of  bulk 
(subscript  B)  and  of  the  grain  boundary  (subscript  GB)  are 
respectively,  kB.  CB,  Kgb.  and  Cgb  [32],  At  the  electrode/electrolyte 
interface,  electrochemical  reactions  may  take  place  which  results 
in  the  faradaic  impedance  ZF.  In  contrast  to  a  usual  electro¬ 
chemical  system,  the  high  frequency  limit  of  the  impedance  does 
not  tend  to  a  resistance  since  the  resistance  of  leads  and  con¬ 
nectors  are  negligibly  small,  therefore  the  high  frequency  limit  is 
practically  zero.  However,  a  self-inductance  due  to  these  compo¬ 
nents  and  a  stray  capacitance  of  instruments  may  no  longer  be 
neglected  in  the  very  high  frequency  domain,  for  instance  above 
1  MHz.  It  is  worth  to  emphasize  that  many  literature  work  neglect 
this  crucial  aspect  leading  to  an  erroneous  estimation  of  the  re¬ 
sistivity.  As  a  whole,  the  electrical  equivalent  circuit  presented  in 
Fig.  2  is  proposed. 

The  impedance  of  this  equivalent  circuit  can  be  expressed  by: 


7  —  ; .  c . , ,  _i _ YE _ 1 _ ub _ 

J  +1+0'^B'CB)  1+a-^GB  CGB) 

+ _ A _ 

Ri+ZJ+  (j'w'ftt'Cd)ad 


where  ZF  is  considered,  as  often  reported,  by  the  series  combination 
of  a  charge  transfer  resistance  Rt  and  of  a  diffusion  impedance  Zd.  Cd 
is  the  double  layer  capacitance. 

In  this  model,  the  capacitances  C  (B  or  GB)  are  assumed  to  be 
frequency  dependent  due  to  the  dielectric  relaxation  phenomenon. 
One  of  the  simplest  models  for  describing  this  phenomenon,  the 
Debye  relationship  with  Cole— Cole  type  frequency  distribution, 
was  adopted  [56], 


L  =  length 


grain 

boundaries 


Electrode 
(Area  =  A) 

G  -  g  =  Grain 

Edge  length  boundary  thickness 


Fig.  1.  Brick  Layer  Model  (BLM)  to  analyse  impedance  data  for  of  solid  electrolyte. 
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Fig.  2.  Equivalent  electrical  circuit  based  on  BLM  for  fitting  the  impedance  results;  S 
represents  a  spurious  inductance  due  to  the  cabling;  the  left  hand  circuit  represents 
the  bulk  (B),  the  center  one  the  grain  boundary  (GB),  and  the  right  hand  one  for  the 
electrode/electrolyte  interface  (current  injection  electrodes). 


C  —  Coo  + 


Co  -  C« 
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The  parameters  stand  for: 


(3) 


Cra:  High  frequency  limit  of  the  capacitance  (B  or  GB)  (F), 
C0:  Low  frequency  limit  of  the  capacitance  (B  or  GB)  (F), 
j:  Imaginary  number;  V—T, 
t :  Relaxation  time  constant  (s), 
a:  Cole  —  Cole  coefficient;  0—1  (dimensionless). 


The  capacitance  C  and  the  relative  permittivity  e  of  the  ceramic 
are  defined  by  the  planar  condenser  model: 


C  = 


e-ey-A 


where, 


(4) 


ey:  Dielectric  constant  of  vacuum;  8.85  10~14  (F  citT1), 
A:  Geometrical  electrode  surface  area  (cm2) 
d :  Thickness  of  the  electrolyte  (cm) 


The  diffusion  impedance  can  be  derived  from  the  Warburg 
expression  for  a  bounded  system  and  expressed  as: 


2d  =  Rd 


tanh^/(j-6j-rd)p__  Rd 

y/CjAAd)p  (l+j-wTd) 


(5) 


The  frequency  dependent  term  of  the  most  right  hand  term  in 
Eq.  (5),  Cole-Davidson  type  dispersion  [57],  is  easier  for  computer 
calculation  than  the  central  expression,  and  many  commercially 
available  softwares  do  not  handle  this  equation.  If  =  1,  then  Zd  is 
practically  identical  to  the  bounded  diffusion  Warburg  impedance. 

Simulation  of  Eq.  (2)  displayed  in  Fig.  3  then  represents  the 
impedance  of  a  proton-conducting  electrolyte  without  the  contri¬ 
bution  of  the  spurious  inductance  S. 

Instead  of  a  depressed  semicircle,  an  asymmetry  can  be  seen  for 
ZB  and  Zgb,  and  this  is  due  to  the  capacitance  distribution  repre¬ 
sented  by  Eq.  (3). 

The  values  of  the  parameters  used  for  the  simulation  are  rep¬ 
resented  in  Table  1. 

With  the  BLM,  the  number  of  elemental  cube  layers  n  can  be 
calculated  by: 


n  = 


L  _A 
g  +  G~G 


(6) 


0  100  200  300  400  500  600 

Real  part  /  Q 


Fig.  3.  Simulated  impedance  spectra  of  a  perovskite  (cf.  Fig.  2);  ZT:  Overall  impedance 
spectrum;  Zb:  Impedance  relative  to  the  bulk;  Zgb:  Impedance  at  grain  boundaries;  Zp: 
Faradic  impedance  due  to  the  electrode  response.  The  contribution  of  self-inductance  S 
by  leads  and  connectors  is  neglected. 


1  £B£VG2  1  r- 
Cb  =  n  G  =  n'£B'tv 


(7) 


Similarly,  the  capacitance  related  to  the  grain  boundary  Cgb  will 
be  evaluated  by: 


Cgb  =  - 


1  eGB-wG2 


(8) 


Provided  that  the  relative  dielectric  constant  at  high  frequency 
limit  of  the  bulk  eB,oo  and  of  the  grain  boundary  egb.oo  are  identical, 
then  it  is  readily  derived  from  Eqs.  (6)  and  (7): 


O) 

'-GB,oo  '-J 

Similarly,  one  can  evaluate  the  conductivity  of  bulk  <rB  and  grain 
boundaries  ogb  as  follows: 


1  L 
~  Rb  A 


(10) 


ffGB 


Rgb  4  G 


(11) 


The  parameter  values  of  the  equivalent  circuit  were  determined 
by  a  house-made  simplex  regression  calculation,  and  then  <rB  and 
ogb  were  determined  from  regressed  data. 


3.  Experimental 

Ln-substituted  SrZro.9Lno.1O2.95  ceramic  powders  with  an 
average  grain  size  of  30  nm  were  prepared  as  reported  in  Refs. 
[49—541.  The  synthesized  powders  were  then  sintered  to  a  pellet 
shape  of  1.5  cm  in  diameter  and  ca.  3  mm  in  thickness  by  one  of  the 
two  sintering  methods: 

-  Conventional  pressing  into  disks  followed  by  a  sintering  at 
1480  °C  under  40  MPa  with  additives  for  2  h  in  an  argon 
atmosphere. 


Table  1 

Values  of  variables  used  for  the  simulation  calculations. 


ZB  Rg  =  150  Q  Cg.oo  =  0.1  nF  CG>0  =  0.16nF 

Zgb  ^gb  =  300  Q  Cgb.oo  =  ^  ^gb.o  =  40  nF 

Zp  Rt  =  20  Q  Q  =  60  pF  ad  =  0.7 

Ri  =  oo  Q  Ci  =  20  pF  /?  =  0.66 


The  capacitance  Cb  can  then  be  evaluated  with  Eq.  (2): 


(Xq  —  1.00 

aGB  —  0.75 


0.  Lacroix  et  al.  /  Journal  of  Power  Sources  270  (2014)  506—515 


509 


-  Uniaxial  hot  pressing  into  disks  between  40  and  1000  MPa  in  air 
at  1600  °C. 

Elementary  analyses  revealed  that,  according  to  the  synthesis 
process  and  ionic  radius  of  the  substituent  element,  the  value  of  x 
varies  between  0  and  0.1  in  the  general  perovskite  formula  given 
above.  We  focussed  essentially  on  the  highest  substituted 
perovskites. 

Electrochemical  measurements  were  carried  out  with  a  poten- 
tiostat  (Solartron  1287)  and  a  frequency  response  analyser  (Solar- 
tron  1255).  For  several  measurements  including  very  high 
frequency  domain  (up  to  10  MHz),  we  used  the  Solartron  1260.  The 
behaviour  of  samples  as  disks  was  investigated  between  400  and 
800  °C  under  various  atmospheres. 

The  results  obtained  by  thermo-gravimetric  analyses  (TGA) 
were  previously  described  by  Sala  et  al.  [48].  High  temperature 
water  uptake  with  different  hydration  methods  was  then  carried 
out  on  these  ceramics  in  various  conditions.  Electrochemical  and 
thermo-gravimetric  measurements  were  then  carried  out  to  these 
samples: 

-  Electrochemical  measurements  were  performed  in  a  furnace 
with  a  sweeping  of  wet  argon  (with  10  vol-%  of  steam)  or  in  an 
autoclave  operating  at  0—10  bars  of  steam  pressure  between 
400  and  800  °C. 

-  Water-uptake  was  estimated  by  thermo-gravimetric  analysis  of 
pre-pressurised  ceramics  in  autoclave.  The  standard  treatment 
of  ceramics  samples  was  200  °C,  steam  pressure  at  15  bars 
(1.5  MPa)  during  36  h.  Other  measurements  were  carried  out 
after  36  h  of  treatment  at  263  °C  and  steam  pressure  at  50  bars 
(5  MPa). 


4.  Results 

The  morphology  of  sintered  perovskite  ceramic  is  presented 
first.  Then,  the  Thermo-Gravimetric  Analysis  (TGA)  measurements 
for  characterizing  the  water  molecules  release  thus  the  removal  of 
proton  from  solid  electrolyte  will  be  estimated  qualitatively.  Finally, 
the  characterization  of  perovskite  electrolyte  by  electrochemical 
impedance  spectroscopy  will  be  shown. 

4.3.  Morphological  characterisation 

Fig.  4  shows  the  SEM  picture  of  the  perovskite  surface  sintered 
by  two  different  methods.  The  grain  size  changes  with  the  sintering 
method.  Conventional  sintering  with  additive  presents  an  average 


Temperature,  T I  °C 

Fig.  5.  Thermo-gravimetric  analyses  of  a  substituted  perovkite.  a)  Standard  hydration; 
d,i,  =  88%;  ATG  under  Ar;  CP  (G  =  1  pm),  b)  Standard  hydration;  dth>98%;  ATG  under 
02;  UP  (G  =  300  nm). 


grain  size  of  1  pm  whereas  by  uniaxial  hot  pressing,  grain  size  is 
close  to  300  nm.  These  two  materials  will  be  labelled  further  as  CP 
(“Conventional  Pressing)  and  UP  (“Uniaxial  Hot  Pressing”).  High 
densification  rates  of  98—100%  were  achieved  by  both  methods. 

4.2.  Gravimetric  measurements 

As  previously  reported  [48]  the  TGA  of  substituted  perovskite 
conditioned  by  the  standard  hydration  treatment,  exhibits  four  loss 
steps  as  the  temperature  increases  (Fig.  5). 

The  threshold  temperature  of  weight  loss  is  determined  by  the 
inflexion  points  that  indicate  different  degrees  of  interactions  be¬ 
tween  the  inserted  water  and  the  ceramic.  The  first  two  weight  loss 
steps  (around  50  and  150  °C)  are  attributed  to  inserted  water  in 
ceramic.  The  third  and  the  fourth  weight  loss  steps,  around 
400—800  °C  are  assigned  to  the  loss  of  thermally  stable  species 
such  as  H+  or  OH~,  depending  on  their  type  of  “hydrogen  bonding” 
with  neighbouring  oxygen. 

The  comparison  of  curves  “a”  and  “b"  indicates  that  the  mass 
loss  is  more  marked  for  the  UP  ceramic.  Other  experiments,  not 
presented  here,  showed  more  noticeable  mass  loss  for  less  densi- 
fied  ceramics.  Above  800  °C,  and  under  oxygen  flow,  a  weight  in¬ 
crease  is  observed  after  water  loss.  Oxygen  species,  such  as  02~,  fill 
the  vacancy  sites  again  (Eq.  (1)).  This  behaviour  confirms  the 


Fig.  4.  SEM  micrographs  of  the  sintered  perovskites;  (a)  Conventional  sintering  (CP)  at  1480  °C  under  40  MPa  with  additives  under  Ar;  (b)  Uniaxial  hot  pressing  (UP)  at  40  MPa  in  air 
at  1600  °C  without  additive. 
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Table  2 

Filling  rate  of  several  pre-pressured  ceramics  calculated  from  thermo-gravimetric 
data. 


Specimens 

Grain 
size  (nm) 

TGA 

atmosphere 

Densification 
rate  (dth%) 

Hydration 

pressure 

Filling 
rate  (%) 

CP 

1000 

n2 

88 

10  bars 

100 

CP 

1000 

o2 

88 

10  bars 

100 

CP 

1000 

02 

>98 

10  bars 

1 

CP 

1000 

02 

>98 

50  bars 

40 

UP 

300 

n2 

>98 

10  bars 

75 

oxygen  mobility.  It  is  worth  noting  that  in  this  case,  ion  conduction 
is  likely  no  longer  controlled  by  protons  but  by  anions  or  vacancies. 

The  degree  of  water-filling  was  evaluated  from  the  ratio  of  the 
relative  mass  loss  to  the  substituent  concentration.  Table  2  sum¬ 
marises  the  filling  rate  together  with  other  experimental  data  for 
different  densification  rates  and  steam  pressures. 

This  table  indicates  that  the  atmosphere  at  which  the  TGA  was 
carried  out  has  no  effect  on  the  measured  water-filling  rate. 

After  a  standard  hydration  treatment  (10  bars),  CP  ceramics 
(grain  size  around  1  pm)  are  completely  filled  up  by  water  for  a 
densification  rate  of  88%  while  only  1%  of  vacancies  are  filled  when 
the  ceramic  densities  are  greater  than  98%.  However,  pressure  in¬ 
crease  overcomes  the  difficulties  of  hydration  on  highly  densified 
ceramic.  An  enhancement  of  steam  pressure  allowed  filling  of  40% 
of  vacancy  sites. 

Under  the  same  hydration  pressure,  for  ceramics  with  the  same 
densification  rate,  it  is  easier  to  fill  the  vacancies  when  the  grain 
size  is  smaller.  As  an  example,  for  the  98%  densification  rate,  the 
filling  rate  rises  from  1  %  to  75%  by  reducing  the  grain  size  from  1  pm 
to  300  nm  (see  Table  2). 

4.3.  Electrochemical  impedance  measurements 

On  both  sides  of  perovskite  pellet  of  1.3  cm  in  diameter  and  ca. 
2  mm  in  thickness,  platinum  meshes  were  inserted  as  electrodes. 
The  impedance  measurements  were  carried  out  at  zero  bias 
current. 

Fig.  6  presents,  as  an  example,  the  impedance  spectrum  ob¬ 
tained  for  the  perovskite  specimen  at  500  °C  under  wet  Ar  (10  vol-% 
steam)  flow.  The  specimen  was  left  in  an  autoclave  at  200  °C  and 
15  bars  during  36  h  to  ensure  the  water  filling. 

In  this  figure,  two  capacitive  loops  are  clearly  seen,  one  related 
to  the  bulk  and  another  to  the  grain  boundary.  These  experimental 
data  validate  the  equivalent  circuit  presented  in  Fig.  2.  However,  in 
most  cases,  the  two  time-constants  are  poorly  distinguished,  and 


Real  Part  /  kn 

Fig.  6.  The  impedance  diagram  of  perovskite  ceramic,  densified  at  88%,  pre-hydrated. 
under  15  bars  at  200  °C  during  36  h,  in  autoclave  and  then  tested  under  wet  argon  flow 
(10  vol-%  Steam)  during  12  h  at  500  °C. 


the  parameter  regression  calculation  is  necessary  to  separate  these 
two  contributions  in  the  overall  conductivity  of  the  ceramics.  It  is 
difficult  to  know,  a  priori,  which  structure,  bulk  or  grain,  appears  at 
higher  frequencies,  since  the  time-constant  is  determined  by  R— C, 
where  R  is  expected  to  be  proportional  to  the  thickness  whereas  the 
capacitance  is  inversely  proportional  to  it.  It  is  expected  however 
that  the  lower  capacitance  obtained  is  due  to  the  bulk,  since  its 
thickness  is  much  larger  than  that  of  the  grain  boundaries.  In  all 
experimental  data  of  this  work,  the  impedance  of  bulk  appeared  in 
the  highest  frequency  domain. 

It  is  important  to  remark  that  the  low  frequency  limit  of  the 
impedance  spectrum  tends  to  infinite  in  an  inert  gas  atmosphere. 
This  is  the  signature  of  the  absence,  as  expected,  of  faradaic  reaction 
but  also  of  electronic  conduction  through  the  ceramic. 

4.3.1.  Reversibility  of  water  insertion 

As  mentioned  above,  since  the  perovskite  ceramics  are  not 
intrinsically  proton  conducting  materials,  the  insertion  of  water  in 
their  structure  is  necessary.  It  is  expected  therefore  that  the  ionic 
conduction  depends  on  the  amount  of  water  having  interacted  with 
the  solid  electrolyte.  For  sake  of  comparison,  the  results  are  pre¬ 
sented  in  a  normalized  scale,  i.e.  plotting  the  R  ratio  between  the 
ionic  conducting  resistance  (Rg  +  Rgb)  at  time  t  and  its  steady-state 
value  at  the  wet  state. 

Rs  =  Rb  +  Rgb  (12) 


r=mi 

K£,wet 


(13) 


Fig.  7  depicts  the  variation  of  the  normalized  resistance  R  for  a 
CP  sintered  ceramic  densified  at  88%. 

The  experiments  were  carried  out  as  follows.  First,  the  ceramic 
specimen  was  held  during  24  h  under  He  environment  at  the 
pressure  of  10  bars  at  a  given  temperature  (500  or  600  °C  for  the 
example  illustrated  here).  Then,  He  gas  was  replaced  by  steam  at  a 
pressure  of  10  bars.  The  impedance  at  1  kHz  was  recorded  every 
30  s  to  follow  its  evolution.  R2  was  evaluated  from  the  impedance 
modulus  and  the  phase  shift.  After  several  hours  under  this  at¬ 
mosphere,  depending  on  the  temperature,  He  gas,  at  the  same 
10  bars  pressure,  was  introduced  into  the  electrolysis  cell. 

For  both  temperatures,  as  soon  as  the  ceramic  is  exposed  to 
steam  atmosphere,  a  steep  decrease  of  the  ionic  resistance  was 
observed,  In  contrast,  when  He  gas  replaces  steam  in  the  electrol¬ 
ysis  cell,  R  increases  slowly,  and  regularly  with  time.  Two  incidents 
are  however  observed  for  this  monotonic  R  increase  (around  4  h  at 
500  °C  and  70  h  at  600  °C).  These  events  are  due  to  a  small  pressure 
drift  during  the  experiment,  which,  once  corrected,  allowed  a 
satisfactory  extrapolation  of  the  initial  part  of  the  curves.  In  this 
figure,  it  can  be  remarked  also  that  only  10  h  are  necessary  to 
remove  the  major  part  of  the  inserted  water  at  500  °C  while 
60  h  are  needed  at  600  °C.  One  could  explain  this  paradoxical  result 
by  a  structural  symmetry  distortion  from  orthorombic  to  cubic 
when  the  temperature  changes  from  500  °C  to  600  °C.  This  sym¬ 
metry  change  induces  the  modification  of  local  interaction  with 
charge  careers  [58,59],  This  behaviour  is  crucial  because  it  can  in¬ 
fluence  the  continuity  of  operation  of  the  electrolysis  cell  in  the 
case  of  temperature  dysfunction. 

The  insertion  of  water  molecules  is  checked  also  by  measuring 
the  relative  dielectric  constant.  The  impedance  measurements 
were  carried  out  between  1  MHz  and  10  mHz  every  30  min.  The 
Cole— Cole  capacitance  C(oj)  is  calculated  according  to  the  following 
equation  after  correcting  the  parasitic  inductance  of  the  high  fre¬ 
quency  domain: 
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Fig.  7.  The  evolution  of  the  ionic  resistance  RE  =  RB  +  Rcb  evaluated  from  impedance  measurements  at  10  Hz  under  the  steam  or  He  environment  under  10  bars,  (a)  at  500  °C  and  (b) 
at  600  “C. 


By  extrapolating  Cu  plot  in  the  complex  plane  towards  the  high 
frequency  limit,  and  applying  Eq.  (4),  the  relative  dielectric  constant 
is  calculated.  Fig.  8  presents  the  evolution  of  e/ewet  in  reduced  scale 
where  ewet  stands  for  the  relative  dielectric  constant  of  the  ceramics 
at  the  steady  state  under  the  water  vapour  atmosphere. 

It  can  be  seen  that  the  dielectric  constant  increases  during  the 
insertion  of  water  molecules  in  the  ceramic  structure,  and  de¬ 
creases  when  the  content  of  water  molecules  diminishes  under  He 
environment. 

Figs.  7  and  8  indicate  that  the  insertion  of  water  molecules  is 
reversible  but  the  output  of  water  molecules  into  the  ceramic 
structure  is  a  rather  slow  process.  The  difference  is  due  to  the  re¬ 
action  of  water  with  oxide  ion  vacancies  of  ceramic  matrix.  The 
protonic  species  are  stabilised  by  the  interaction  with  dopant 
element,  which  have  an  oxidation  degree  lower  than  zirconium 
site.  Thus,  local  interaction  is  modified  due  to  dopant,  then  the 
removal  of  water  becomes  more  difficult  than  insertion  process 
[58,59],  Besides,  since  the  anodic  reaction  consumes  the  water 
molecules,  the  electrolyte  should  be  kept  wet  to  maintain  a  high 
ionic  conduction. 

4.3.2.  Effect  of  temperature  on  relative  permittivity 

Fig.  9  illustrates  the  variation  of  permittivity  with  temperature 
calculated  from  Cb.oo  with  Eqs.  (8)  and  (2). 

With  perovskite  ceramic,  densified  at  95%,  the  relative  permit¬ 
tivity  increases  with  temperature,  and  also  under  steam  flow.  At 


Fig.  8.  Effect  of  inserted  water  molecules  on  the  dielectric  constant  at  500  and  600  °C. 
Steam  or  He  pressure  at  10  bars;  reduced  value  with  respect  to  the  wet  sate. 


600  °C  for  instance,  is  close  to  100  under  air  atmosphere  and  to 
350  under  steam  flow.  This  high  permittivity  would  reflect  the 
modification  of  material  due  to  water  dissociation  and  formation  of 
hydroxyl  ion  in  the  material  as  seen  in  Eq.  (1). 

As  seen  above,  the  evolution  of  permittivity  from  a  dry  state  to  a 
hydrated  state  seems  to  be  linked  to  hydroxyl  ion  formation.  Thus, 
the  comparison  of  CBi00  and  Cgb.o,  can  also  give  an  idea  about  the 
location  of  the  steam  insertion  (Table  3).  As  explained  before,  under 
dry  conditions,  Eq.  (9)  gives  an  estimate  of  the  geometrical  pa¬ 
rameters  of  the  ceramic  structure. 

As  the  Cb.oo/Cgb.oo  ratio  is  proportional  to  the  grain  boundary 
thickness  (g/G),  this  could  be  explained  by  the  fact  that  the  use  of  a 
sintering  additive  thickens  the  grain  boundary. 

The  evolution  of  this  ratio  under  steam  treatment  indicates 
opposite  behaviours  for  the  two  ceramics.  For  the  UP  ceramic  under 
steam,  the  ratio  increases  while  it  decreases  for  the  CP  ceramic.  It 
can  be  concluded  that  water  insertion  depends  on  grain  boundary 
size,  on  the  distribution  of  vacancies  and/or  on  the  sintering  ad¬ 
ditive  content. 

4.3.3.  Effect  of  steam  pressure 

It  was  observed  that  the  insertion  of  water  molecules  into  the 
ceramic  structure  plays  a  determining  role  on  the  ionic  conduc¬ 
tivity.  Fig.  10  presents  the  evolution  of  the  impedance  spectrum  in 
absence  of  water  and  under  different  steam  pressures  at  600  °C 

[60], 


Fig.  9.  Evolution  of  the  dielectric  constant  estimated  from  the  high  frequency  limit  of 
the  impedance  as  calculated  from  Eq.  (4)  in  two  atmospheres  for  a  ceramic  with  95%  of 
densification  rate. 
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Table  3 

Mean  Cgb/Cb  ratio  between  400  and  700  °C  in  dry  and  wet  conditions  for  same 
powder  sintered  under  various  conditions. 


Sintering  (densification  rate) 

Grain 

Cgb/Cb 

Cgb/Cb 

size  (nm) 

under  Ar 

under 

steam 

Conventional  with  additive  (>98%) 

1000 

9 

5 

Uniaxial  Hot  Sintering  (>98%) 

300 

15.7 

62 

In  this  figure,  it  can  be  noticed  that  the  overall  ionic  conduction 
increases  with  the  steam  pressure,  which  can  be  explained  by  an 
increase  of  the  water  insertion  in  the  ceramic  and  therefore  of  the 
number  of  charged  carriers  generated. 

In  Fig.  11  the  overall  conductivity  is  displayed  with  respect  to 
temperature  measured  at  different  steam  pressures.  It  can  be  seen, 
as  expected  that  the  higher  is  the  temperature,  and  the  steam 
pressure,  the  higher  is  the  overall  conductivity  of  the  solid 
electrolyte. 

As  an  example,  for  a  steam  pressure  of  10  bars,  the  overall 
conductivity  is  1  mS  cur1  at  600  °C,  the  highest  temperature 
compatible  with  a  usual  structural  material. 

An  inductive  behaviour  can  be  frequently  observed  at  high  fre¬ 
quencies,  as  shown  in  Fig.  10.  This  spurious  additional  element  is 
due  to  the  cell  configuration  including  wiring  and  electrodes,  and  to 
the  frequency  response  analyser  used.  It  is  worth  to  notice  that  the 
use  of  high  temperature  cell  needs  long  electrical  leads  between 
the  ceramic  specimen  and  the  frequency  response  analyser,  thus 
leading  to  a  rather  high  inductance.  The  parameter  regression 
calculation  was  really  necessary  to  remove  this  parasitic  inductive 
response,  so  as  to  separate  and  correctly  assess  the  contributions  of 
the  bulk  and  the  grain  boundary  in  the  ionic  conduction,  according 
to  the  equivalent  circuit  presented  in  Fig.  2. 

By  using  the  G/g  ratio  evaluated  from  Eq.  (9)  in  dry  conditions, 
the  conductivity  of  grain  and  grain  boundaries  were  evaluated.  The 
results  are  reported  in  Fig.  12  for  two  ceramic  specimens  of 
different  grain  size  (CP  and  UP  ceramics). 

It  can  be  seen  that  the  (bulk  conductivity)  is  much  higher 
than  that  of  grain-boundaries  (<tgb)-  The  bulk  conductivity  of  CP  and 
UP  ceramics  are  similar  especially  for  temperatures  higher  than 
550  °C  which  sounds  consistent.  In  contrast,  the  UP  ceramic  ex¬ 
hibits  higher  grain-boundary  conductivity. 


Fig.  10.  Impedance  spectra  of  UP  ceramic  densified  at  98%  under  argon  atmosphere  at 
10  bars  and  different  steam  pressures  at  600  °C. 


Temperature,  T /  °C 

Fig.  11.  Overall  conductivity  for  UP  under  different  steam  pressure. 

4.3.4.  Activation  energy  of  ionic  conduction 

Temperature  influence  on  conductivity  can  be  displayed  on  an 
Arrhenius  plot. 

a-T  =  Aexp(j^j  (15) 

Fig.  13  presents  the  results  for  bulk  and  grain  boundaries  con¬ 
duction  at  the  steam  pressure  of  10  bars  for  two  ceramics.  It  can  be 
seen  that  the  bulk  conduction  exhibits  three  domains  in  agreement 
with  thermo-gravimetric  results.  In  contrast,  the  grain  boundary 
conduction  does  not  show  clearly  the  intermediate  temperature 
domain. 

The  apparent  activation  energy  calculated  from  these  data  is 
summarized  in  Table  4  for  the  different  temperature  domains  for 
bulk  and  grain  boundaries. 

Except  for  the  intermediate  temperature  domain,  £B  of 
0.44—0.93  eV,  falls  in  a  reasonable  range  for  an  ionic  conductivity.  It 
can  be  observed  also  that  the  activation  energies  at  the  grain 
boundaries  in  the  low  temperature  range  are  higher  than  those 
found  for  the  bulk.  Therefore,  there  is  some  difference  in  the  nature 
of  the  ionic  conductivities  of  the  bulk  and  the  grain  boundary.  It  is 
also  remarked  that  at  the  intermediate  temperature,  the  activation 
energy  is  very  small,  i.e.  the  conductivity  change  little.  In  this 


Fig.  12.  Conductivity  of  bulk  and  grain  boundary  for  the  conventional  sintering  (CP) 
and  uniaxial  hot  pressing  (UP)  ceramics  under  10  bars  of  steam  pressure. 
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Fig.  13.  Arrhenius  plot  of  bulk  and  grain  boundary  conduction  for  conventional 
pressed  (CP)  and  uniaxial  hot  pressed  (UP)  perovskitetype  ceramics  at  steam  pressure 
of  10  bars. 

temperature  range,  however  the  water  molecules  are  removed 
from  the  ceramic  matrix,  thus  the  density  of  charge  careers  de¬ 
creases,  and  also  it  is  very  likely  that  the  ionic  careers  change  from 
protons  to  anions.  Therefore,  the  calculated  activation  energy 
drawn  from  these  data  is  only  apparent.  Except  this  particular 
transition  domain,  the  influence  of  the  ceramic  structure  can  be 
excluded  since  the  conduction  of  the  CP  and  UP  ceramics  are 
essentially  the  same. 

5.  Discussion 

From  the  results  presented  above,  it  can  be  concluded  first  that 
densification  is  one  of  the  most  crucial  parameters  for  optimization. 
High  densification  of  the  material  reduces  the  hydration  of  ce¬ 
ramics  (Table  2;  Fig.  5).  Conversely  the  presence  of  porosity  facili¬ 
tates  the  diffusion  of  water  inside  the  material,  and  also  the  dead¬ 
end  pores  inside  the  ceramic  matrix  may  create  host  sites  for 
inserted  water  molecules.  The  differences  in  hydration  behaviours 
observed  in  ceramics  with  different  grain  sizes  are  a  result  of  the 
grain  boundaries  length  in  the  ceramics.  In  other  words,  the  grain 
boundaries  in  the  material  serve  as  pathways  for  water  distribution 
into  the  ceramic  structure. 

If  it  is  assumed  that  the  diffusion  of  water  is  fast  through  the 
grain  boundaries,  then,  the  rate-determining  step  of  full  water 
insertion  is  essentially  determined  by  the  grain  size.  Therefore,  the 
time-constant  for  hydration  of  the  grain  size  G  and  the  diffusion 
coefficient  of  water  Du.l()  can  be  estimated  by: 


Table  4 

Activation  energy  of  bulk  £B  and  grain  boundary  EGb  between  400  and  800  °C  at 
10  bars  of  steam  pressure  for  conventional  sintered  and  uniaxial  hot  pressed 
ceramic. 


T.  Domain 

T.  Range  (B) 

Eb  (eV) 

Egb  (eV) 

Conventional  sintered 

Low 

400  °C— 550  °C 

0.62 

0.8 

Intermediate 

550  °C— 650 °C 

0.09 

- 

High 

650  °C  - 

0.93 

- 

Uniaxial  hot  pressed 

Low 

400  °C— 550  °C 

0.63 

0.8 

Intermediate 

550  °C— 650 °C 

0.08 

- 

High 

650  °C  - 

0.44 

0.45 

This  equation  suggests  that  the  penetration  of  proton  defects 
could  be  accelerated  by  an  increase  of  water  pressure  (Table  1). 
Besides,  higher  pressure  increases  the  flux  of  water  due  to  the  in¬ 
verse  osmosis,  and  thus,  the  number  of  charge  carriers.  Although 
water  diffusivity  increases  with  temperature  in  the  grain,  the  water 
depletion  is  taking  place  simultaneously  as  attested  by  Fig.  5.  This 
explains  a  small  change  of  the  ionic  conduction  in  the  bulk.  In 
contrast,  no  such  effect  was  observed  at  the  grain  boundaries.  It 
was  shown  that  the  release  of  water  molecules  from  the  ceramics  is 
faster  at  the  bulk  than  at  the  grain  boundaries. 

Fig.  5  shows  that  at  higher  temperatures  (>750  °C),  the  amount 
of  absorbed  water  tends  towards  zero,  whereas  the  ionic  conduc¬ 
tivity  increases.  This  apparent  contradiction  can  by  explained  by 
assuming  a  change  in  the  conduction  nature;  the  protonic  charge 
carrier  being  replaced  by  02~  [61,62],  It  is  worth  reminding  that 
perovskite  is  an  intrinsically  anion  conductor  in  presence  of 
vacancies. 

As  previously  established,  the  conductivity  of  the  bulk  is  higher 
than  that  of  grain  boundary  (see  Fig.  12).  However,  it  is  noteworthy 
that,  even  if  the  conductivity  is  lower  at  grain  boundaries,  its 
contribution  (1/J^gb)  in  the  overall  conduction  (1  /Rs)  is  similar  to 
that  of  bulk  (1/RB)-  This  is  due  to  the  fact  that  the  grain  boundaries 
are  much  thinner  than  the  bulk. 

The  reason  for  a  lower  ionic  conductivity  at  grain  boundaries 
may  be  explained  as  follows:  charged  species  (ions  of  perovskites 
and  impurities),  which  tend  to  segregate  essentially  near  the  grain 
boundaries  neutralize  the  vacancy  sites  at  their  vicinity  for 
ensuring  electro-neutrality.  [63,64],  Consequently,  grain  boundary 
conductivity  decreases  as  reported  for  zirconia  [63  .  It  must  be  also 
emphasized  that  grain  boundary  conductivity  depends  on  sintering 
conditions,  while  no  significant  difference  is  observed  in  bulk 
conduction  (Fig.  13).  Experimentally,  the  ratio  Cgb/Cg  increases  for 
UP  under  steam,  while  it  decreases  for  CP  (Table  3).  There  are  two 
alternative  explanations  about  this  fact.  First,  good  densification 
ratio  was  achieved  by  using  an  additive.  This  element  segregates 
preferentially  at  the  grain  boundary  region  during  sintering  [38] 
leading  to  a  local  decrease  of  ionic  conductivity,  as  this  is  the  case 
for  silica  in  zirconia  [39],  Segregated  species,  such  as  impurities  and 
additives,  are  more  diluted  in  UP  than  in  CP  because  of  a  greater 
grain  boundary  area.  Then,  precipitation  of  elements  at  the  grain 
boundaries,  which  induces,  as  mentioned  above,  a  reduction  of 
vacancy  concentration,  is  smaller  for  UP.  The  effect  of  segregation 
becomes  therefore  much  less  marked  for  uniaxial  hot  pressed 
material. 

Consequently,  there  are  two  possible  ways  to  optimizing  the 
proton  conduction  [65],  First,  since  the  grain  conduction  is  higher 
than  that  of  grain  boundaries,  the  increase  in  grain  size  favours  the 
overall  conductivity.  Ideally,  a  fully  hydrated  single  crystal  might 
exhibit  the  best  conductance.  However,  as  segregations  precipi¬ 
tated  near  the  grain  boundaries  decrease  the  ionic  conduction,  the 
ceramic  should  be  prepared  with  raw  materials  containing  a  low 
concentration  of  impurities  that  decrease  the  vacancy  concentra¬ 
tions.  Second,  it  is  observed  also  that  grain  boundaries  act  as  a 
preferential  pathway  for  water  insertion.  A  smaller  grain  size  fa¬ 
cilitates  therefore  the  insertion  of  water  in  the  bulk  of  ceramics,  and 
decrease  the  local  impurity  concentration.  It  results  that  an  opti¬ 
mum  should  be  found  by  adjusting  these  two  opposite  effects  and 
also  raw  materials  having  low  concentration  of  harmful  impurities, 
such  as  silica,  must  be  selected. 

6.  Conclusion 

The  main  results  obtained  in  the  present  work  are  the  following: 

Water  insertion  is  reversible  in  perovskite  structure  but  it  is  a 
slow  process,  which  needs  several  hours  to  reach  the  steady  state. 
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The  filling  ratio  depends  besides  on  the  steam  pressure  and  on  the 
densification  factor.  The  higher  densification  is,  the  more  difficult  is 
the  filling  at  the  vacancy  sites  by  water  molecules. 

The  filling  ratio  has  a  direct  impact  on  material  conductivity: 
conductivity  increases  with  the  amount  of  protonic  species  in  the 
ceramics  and  with  their  mobility.  Steam  pressure  and  temperature 
have  an  opposite  effect  on  the  proton  species  amount.  The  steam 
pressure  increases  the  population  of  charge  carriers  whereas 
temperature  decreases  it  by  drying  the  ceramic  electrolyte. 
Simultaneously,  increasing  pressure  and  temperature  favours  the 
mobility  of  charged  carriers.  As  a  whole,  conductivity  of  the  bulk 
remains  almost  constant  between  550  and  650  °C. 

The  bulk  conductivity  is  higher  than  the  grain  boundary  one; 
therefore,  a  large  grain  size  may  favour  the  ionic  conduction  of 
perovskite  ceramic.  However,  a  large  grain  boundary  area  enhances 
the  insertion  of  water  molecules  into  the  solid  electrolyte,  because 
the  grain  boundaries  constitute  pathways  for  the  water  molecules. 
Besides,  the  precipitation  of  impurities  that  decrease  the  conduc¬ 
tivity  will  be  lowered  for  a  large  grain  boundary  area.  Therefore,  a 
small  grain  size  increases  the  overall  conductivity.  In  order  to 
optimise  solid  electrolyte  conductivity  it  is  necessary  to  decrease 
the  impurities  of  raw  chemicals  from  which  the  sintering  process 
will  be  performed,  and  adjust  the  grain  size. 

The  steam  pressure  is  one  of  the  most  important  factors  to 
improve  the  proton  conduction  through  the  perovskite  ceramics,  in 
our  experimental  conditions;  the  highest  conductivity  was  ob¬ 
tained  at  the  highest  steam  pressure  examined,  that  is,  at  10  bars. 
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List  of  symbols 

[  ]  concentration  (mole  crrT3) 

A  sample  area  (cm2) 

C  capacitance  (F) 

Co  capacitance  at  low  frequency  limit  for  Debye-type 

distribution  (F) 

Coo  capacitance  at  high  frequency  limit  for  Debye-type 
distribution  (F) 

DH2o  diffusion  coefficient  of  water  (cm2  s~3) 

£a  activation  energy  (eV  moH1) 

/  frequency  (Hz) 

g  grain  boundary  thickness  (cm) 

G  grain  edge  length  (cm) 

L  sample  length  (cm) 

L  inductance  (H) 

Ph2o  pressure  of  steam  (bars) 

Rb  ionic  conduction  resistance  through  bulk  (fi) 

Rgb  ionic  conduction  resistance  through  grain  boundary  (O) 
Rc  charge  transfer  resistance  (O) 

Z  impedance  (fi) 

Zd  diffusion  impedance  (0) 

an  Cole  Cole  distribution  factor  for  circuit  Rn/Cn 

/}  Cole  Davidson  distribution  factor 

£V  permittivity  of  vacuum 

e  relative  permittivity  of  ceramic  (Eq.  (4)) 


e wet  relative  permittivity  of  ceramic  in  the  wet  state 

<tb  bulk  conductivity  (S  cm-1) 

®gb  grain  boundary  conductivity  (S  cm-1) 

r  time  constant  (s) 

to  pulsation  (radian  s_1  =  2  it  f) 
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